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ABSTRACT 

We present ground-based high resolution N-band spectra (At) = lSkms" 1 ) of pure rotational lines 
of water vapor in two protoplanetary disks surrounding the pre- main sequence stars AS 205N and RNO 
90, selected based on detections of rotational water lines by the Spitzer IRS. Using VISIR on the Very 
Large Telescope, we spectrally resolve individual lines and show that they have widths of 30-60 km s -1 , 
consistent with an origin in Keplerian disks at radii of ~ 1 AU. The water lines have similar widths 
to those of the CO at 4.67 /jm, indicating that the mid-infrared water lines trace similar radii. The 
rotational temperatures of the water are 540 and 600 K in the two disks, respectively. However, the 
lines ratios show evidence of non-LTE excitation, with low-excitation line fluxes being over-predicted 
by 2-dimensional disk LTE models. Due to the limited number of observed lines and the non-LTE 
line ratios, an accurate measure of the water ortho/para ratio is not available, but a best estimate for 
AS 205N is ortho/para = 4.5 ± 1.0, apparently ruling out a low-temperature origin of the water. The 
spectra demonstrate that high resolution spectroscopy of rotational water lines is feasible from the 
ground, and further that ground-based high resolution spectroscopy is likely to significantly improve 
our understanding of the inner disk chemistry recently revealed by recent Spitzer observations. 
Subject headings: astrochemistry — planetary systems: protoplanetary disks 



1. INTRODUCTION 

The study of the dynamics and chemistry of the inner 
(R < 10 AU) regions of disks around young pre-main se- 
quence stars is a field in rapid growth, to a large extent 
thanks to the availability of sensitive infrared instrumen- 
tation, including the InfraRed Spectrograph (IRS) on the 
Spitzer Space Telescope. The molecular gas at these radii 
plays a key role in the process of planet formation, not 
only as a necessary reservoir for the formation of gas gi- 
ants, but also in the generation and eventual delivery of 
volatile molecular species to terrestrial planets. 

Recent observations have demonstrated that the mid- 
infrared wavelength range in protoplanetary disks con- 
tains a forest of molecular emission lines, many of them 
due to pure rotational transitions of water vapor (Carr 
& Najita 2008; Salyk et al. 2008; Pascucci et al. 2009; 
Pontoppidan et al. 2010). While Carr et al. (2004) de- 
tected and spectrally resolved rovibrational water lines at 
2.3 /im and Salyk et al. (2008) at 3 /im from protoplan- 
etary disks, the observations of the mid-infrared pure 
rotational lines were all carried out at R = A/AA < 600, 
resulting in unresolved and blended line spectra. The 
rovibrational 2-3 fim lines have much higher excitation 
temperatures, and may not probe the same gas as the 
mid-infrared rotational lines. Further, they also suffer 
from significant line blending, making it difficult to study 
the details of the excitation structure. In this Letter, we 
demonstrate that high resolution spectroscopy of pure ro- 
tational transitions of water in the atmospheric V-band 
window is possible. The rotational water lines observable 
from the ground have excitation temperatures of 3000- 
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6000 K, directly comparable to those of the well-studied 
CO rovibrational lines in the 4.7 /im fundamental band 
(Najita et al. 2003; Blake & Boogert 2004; Brittain et al. 
2007; Pontoppidan et al. 2008; Salyk et al. 2009), allow- 
ing for more accurate relative abundance measurements. 
Further, the lines are well separated from one another 
and do not suffer from blending with other species. 

In this Letter, we present high resolution (R ~ 20 000) 
V-band spectra of AS 205N and RNO 90, disks known 
to have strong water vapor emission (Salyk et al. 2008; 
Pontoppidan et al. 2010), obtained using the Very Large 
Telescope Imager and Spectrometer for the mid-InfraRed 
(VISIR, Lagage et al. 2004). Both disks are likely part of 
the Ophiuchus star forming cloud at a distance of 120 pc 
(Loinard et al. 2008). With stellar luminosities of 4.0 and 
3.5 Lq (Chen et al. 1995) and spectral types of K5 and 
G5, the central stars have likely masses close to 1 M@ 
(e.g., Siess et al. 1997; Salyk et al. 2008; Pontoppidan & 
Blake 2010). RNO 90 in particular, being a single star, 
may be a fairly close analog to the young Sun, while AS 
205N is the primary of a l'/3 binary (160 AU). 

2. OBSERVATIONS 

Six water lines with excitation energies varying from 
from 3300-5800 K were selected using the list of Pontop- 
pidan et al. (2009) in combination with the limitations 
imposed by the availability of filters on VISIR. While 
telluric water generates highly pressure-broadened ab- 
sorption, we find that even from a relatively low site, 
such as Paranal Observatory (2635 m), the atmospheric 
transmission for all our targeted lines is >30%, and often 
significantly better. 

Spectra centered on 12.405 and 12.454 //m were ob- 
served using the cross-dispersed mode of VISIR, with a 
chop-throw of 8". In this mode, one source position is off 
the slit because the minimum chop-throw is larger than 
the slit length. Additional settings centered on 12.829 
and 12.893 /im were observed using the long-slit mode, 
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Fig. 1. — VISIR spectra of rotational water lines observed in AS 
205N. The curves below the spectra show the standard star spec- 
tra, which are a combination of the spectral response function and 
the atmospheric transmission. The velocity range of the spectra 
is referenced to the local standard of rest. The red curve is the 
CO rovibrational (Ad = 1) line shape (J < 8), as observed with 
CRIRES (Pontoppidan & Blake 2010), whose center defines the 
vertical dotted line. The blue curves are lines calculated using an 
LTE 2D disk model (see text for details). 



in which all source positions are on slit, doubling the ef- 
fective exposure time. The effective on source exposure 
times were 1250 s and 2500 s for the cross-dispersed and 
long-slit modes for AS 205N, respectively. For RNO 90, 
these exposure times were doubled. All observations were 
carried out using the 0'.'75 slit during good seeing con- 
ditions, specifically ensuring that components of the AS 
205 binary are well separated. AS 205N was observed 
on August 5 and 6, 2009, while RNO 90 was observed 
on separate nights between August 9, 2009 and Septem- 
ber 11, 2009. The spectral images were co- added using 
the ESO VISIR pipeline version 3.2.2. Chopping paral- 
lel to the slit removes most of the sky background, but 
a significant residual, especially from rapidly varying at- 
mospheric water lines, is still visible. These residuals 
were removed by subtracting the median of each row (in 
the cross-dispersion direction) from the spectral image. 
The ID spectra were generated using optimal extraction 
(Home 1986). The wavelength calibration from the ESO 
pipeline was used, and is computed by cross correlating 
the sky emission spectrum with an atmospheric model 
generated using HITRAN (Rothman et al. 2005). 

The spectral response function and atmospheric trans- 
mission spectrum were measured using observations of a 
bright infrared excess source known to not display strong 
water emission at 12 (an in the survey of Pontoppidan 
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Fig. 2. — Same as Figure 1, but for RNO 90. The lower right 
panel presents the average of the three detected lines. The asym- 
metry visible as an increase in flux on the blue side of the detected 
lines is possibly real, but since it coincides with the telluric wa- 
ter lines, we are hesitant to discuss it in detail without further 
observations, preferably with different reflex motion doppler shifts. 



et al. (2010), specifically the bright Hcrbig Ac/Be star 
HD150193. Ideally, a hot star without excess should be 
used, but no such source is sufficiently bright in the N- 
band. One could use a later type star, but in that case 
photospheric absorption lines will be present that may 
be difficult to remove. Hence, the choice was made to 
use an excess source dominated by continuum emission 
and with no apparent lines. Corrections for small shifts 
(of order a few pixels) between the source spectrum and 
that of the standard as well as a correction for differences 
in water column densities using a simple Lambert-Beer 
law were made to minimize the telluric residuals. The 
absolute fluxes were scaled to match those of the Spitzer 
spectra in Pontoppidan et al. (2010). 

For comparison with the water lines, we use archival 
M-band high resolution (R=100,000) spectra of the CO 
rovibrational fundamental (Av = 1) lines observed with 
CRIRES. The CO lines have roughly the same upper 
level energies as the N-band rotational water lines. The 
CRIRES data set will be described in greater detail in a 
spectro-astrometric survey (Pontoppidan & Blake 2010). 
It was processed following the procedures of Pontopp- 
idan et al. (2008). Finally, we extracted an isolated 
lower energy line for comparison from the Pontoppidan 
et al. (2010) Spitzer spectra. The low-resolution (R=600) 
Spitzer spectra do not allow for the measurement of many 
unblended lines, but there are a few exceptions, one of 
them being the IO4 7 — > 9^ e ortho line at 33.5 fim. 
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TABLE 1 

H2O LINE FLUXES 





Line 


Wavel. 


^upper 


AS 205N a 


RNO 90 a 






/im 


K 






17 4 


13 ~> I63 14 (0) 


12.396 


5746 


4.8 ± 0.2 


0.7 ±0.08 


163 


13 -» 152 14 (p) 


12.407 


4915 


2.3 ± 0.2 


0.9 ± 0.11 


13 7 


6 -> 12 4 9 (0) 


12.453 


4187 


10.4 ±0.2 


1.0 ±0.14 


Us 


3 -> 10 5 6 (0) 


12.445 


3606 


11.4 ±0.3 


< 0.3 


125 


7 -» 11a 10 (p) 


12.893 


3290 


1.9 ± 0.2 


< 0.3 


10 8 


2 ~> 95 5 (P) 


12.832 


3223 


< 0.4 


< 0.2 


10 4 


7 -> 93 6 (0) 


33.501 


2260 


1.9 ± 0.2 


< 0.3 



a 10 -14 erg cm -2 s _1 
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Fig. 3.— Rotation diagram for AS 205N and RNO 90. The star 
symbols are ortho lines, while the square symbols indicate para 
lines. The points include a spin degeneracy of 3:1. The dashed 
lines are the best fit single temperature models (using the detected 
VISIR lines only). The units are cgs and the symbols have their 
usual meaning. 

3. RESULTS 

The reduced VISIR spectra of AS 205N and RNO 90 
are shown in Figures 1 and 2, respectively. For AS 205N, 
we detect 5 out of the 6 targeted lines, while 3 lines are 
detected for RNO 90. Integrated line fluxes are given in 
Table 1. 

The lines, where detected, are spectrally resolved, with 
line widths of -SOkms" 1 for AS 205N and -^kms" 1 
for RNO 90. In the figures, the water lines are compared 
to the CO Ai>=l line shape. While variations are appar- 
ent for AS 205N, these differences are relatively minor, 
with FWHM of 20-35 km s" 1 . similar to, or slightly less 
than those of CO. This indicates that the water lines are 
formed at similar or slightly larger disk radii. The high 
energy lines (17 4 3 -> 16 3 4 , I63 13 -> 15 2 14) of AS 205N 
show indications of being flat-topped, while the lower ex- 
citation lines are single peaked, as is CO. All three de- 
tected RNO 90 lines are consistent with being double 
peaked, also matching the corresponding CO lines. The 
double peak is apparent when the three detected line pro- 
files are averaged. We interpret this as evidence for a line 
origin in a Keplerian disk. 

In order to estimate the properties of the emitting 
gas (rotational temperatures, radial distribution, column 
densities and abundances), we model the spectra using 
the 2D raytracer RADLite (Pontoppidan et al. 2009). 
We also construct rotation diagrams (Figure 3), which 
are designed to allow for easy comparison with a model 
consisting of a single-temperature slab of gas, and en- 
able a discussion of the excitation and column density, 
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Fig. 4. — Continuum-subtracted surface brightness profiles of 
three lines for the RNO 90 model. 

in an average sense, of the observed species; LTE level 
populations and optically thin emission predicts that line 
fluxes arrange along a straight line in the rotation dia- 
gram. Here we demonstrate that there are departures 
from the basic slab model assumptions of LTE, and that 
the assumption that all lines are formed in the same 
solid angle does not hold for water. Specifically, the 
slab model predicts significantly brighter 10s 2 ~^ 9s 5 
and 12 5 7 — > H2 10 emission, by at least factors of 5-10, 
than what is observed for AS 205N. The upper limits 
on the RNO 90 lines are not strongly constraining, and 
the rotation diagram for this disk is roughly consistent 
with the slab model in the N-band, while the Spitzer 
33.5 L«m line is undcrluminous. The lower energy Spitzer 
line thus confirms the N-band turnover at lower ener- 
gies in the rotation diagram for AS 205N. We stress that 
these are differences that are not apparent from the low- 
resolution Spitzer spectra alone. A linear regression of 
the VISIR lines on the linear part of the rotation dia- 
gram (Supper > 3500 K) yields rotational temperatures 
of 540 and 600 K for AS 205N and RNO 90, respectively, 
for optically thin emission. 

The AS 205N temperature can be compared to the 
1000 K required to fit the 3 /xm rovibrational lines (Salyk 
et al. 2008), indicating that the rotational lines indeed 
trace cooler gas. The selection of lines includes ortho- 
and para-transitions. The assumption that the N-band 
lines should form a straight line in the rotation diagram 
allows an estimate to be made of the ortho/para ra- 
tio. The paucity of lines, and the presence of clear de- 
partures from LTE and/or slab geometry unfortunately 
makes such an estimate uncertain. Nevertheless, scaling 
the three ortho-lines for AS 205N with energies between 
3500 and 6000 K to an optimal alignment with the two 
detected para lines results in an apparent O/P ratio of 
4.5 ± 1.0 for AS 205N The ratio is higher than 3 because 
the para lines are placed below a straight line defined 
by the three ortho lines, but by less than 2a. While 
this measurement is certainly affected by the caveats dis- 
cussed above, we interpret this as indicative of a warm 
gas-phase (equilibrium) chemistry, as opposed to water 
formed in a low-temperature environment (<60K) via 
grain-surface chemistry, which would produce O/P ratios 
<3. Indeed, cometary water has O/P ratios correspond- 
ing to spin temperatures of ~30K (Crovisier et al. 1997; 
Kawakita et al. 2004), and optical depth effects would 
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drive the estimated 0/P values lower. 

Since the lines are now spectrally resolved, the emitting 
area and hence the column density of the emitting gas 
can be directly estimated by converting the line width 
to a physical location in a Keplerian disk. The incli- 
nation of the RNO 90 disk is -45° (Pontoppidan & 
Blake 2010) and the star has a mass of ~O.9M (An- 
drews et al. 2009). The RNO 90 water lines extend to 
~ 30kms _1 , corresponding to an inner water line emit- 
ting radius of 0.45 AU. Similarly, for AS 205N the line 
width at zero flux corresponds to an inner emitting radius 
of 0.1-0.4 AU, for the 11 8 3 ->• 10 5 6 and 17 4 13 -> I63 14 
lines, respectively, assuming Keplerian flow, a disk in- 
clination of 25° and a stellar mass of 1 Mq (Andrews 
et al. 2009). The effective outer radii of the emitting ar- 
eas are determined by generating 2D models for the lines 
with RADLite, assuming LTE level populations. The 
model uses a simple parametrized flared disk structure 
(Dullcmond et al. 2001), with inner radii fixed accord- 
ing to the observed line widths. The water abundance 
is set to 2.6 x 10~ 4 per H and the disk surface gas-to- 
dust ratio set to 1.28 x 10 4 , simulating significant dust 
settling, as discussed in Meijerink et al. (2009). Finally, 
the stellar luminosities were adjusted to match the dust 
continuum at 12.5 fim. The resulting line strengths and 
profiles match well to those of the detected N-band lines 
with E >4000K, as was already suggested by the rota- 
tion diagram (see Figures 1 and 2). Figure 4 shows the 
corresponding line surface brightness from the RADLite 
model, demonstrating that the RNO 90 emission is dom- 
inated by radii from 0.4 to 1.0 AU (-2.6 AU 2 ), but with 
the emitting area changing by a factor 2-3, depending on 
the line. Due to the higher luminosity of AS 205N, emit- 
ting radii span a wider range, here from 0.35 to 2.5 AU 
(—19 AU 2 ). Using these measures, the column densities 
of water, averaged over the emitting radii, are 2 x 10 18 
and 3 x 10 18 cm~ 2 . 

4. DISCUSSION 

We have detected pure rotational lines from water va- 
por (H2 16 0) in protoplanetary disks around T Tauri 
stars from the ground using transitions in the atmo- 
spheric A-band window. Gas temperatures of 500-600 K 
are found, while the line widths of 30-60 km s -1 confirm 
that the disk surface at radii of 0.4-2.5 AU is forming 
the lines. The resolved lines enable a much more direct 
determination of the emitting area and hence the col- 
umn density of the water, as compared to that possible 
for spectrally unresolved Spitzer observations. We find 
that the assumptions used for the analysis of the Spitzer 
data generally hold, i.e. that the lines are formed in the 
surface of Keplerian disks at radii of ~1 AU. The wa- 
ter column densities are found to be a little higher for 
AS 205N than those determined by Salyk et al. (2008) 
using rotational lines at somewhat longer wavelengths 



and lower energies. This appears to be consistent with 
the rotation diagram that indicates that lower excitation 
lines are less bright than expected for LTE level pop- 
ulations. This effect was also noted by Meijerink et al. 
(2009), who interpreted it as a depiction of water beyond 
~1AU. More extensive N-band studies will help to fur- 
ther constrain the shape of the water rotation diagram to 
constrain the excitation and spatial distribution of water 
vapor in the inner regions of protoplanetary disks. 

The derived water abundance is close to that predicted 
by Bethell & Bergin (2009), based on a chemical model 
that includes water self-shielding. The 2D RADLite 
models of both disks require high water abundances and 
high gas-to-dust ratios in order to reproduce the observed 
lines - a significant fraction of the available oxygen must 
be locked up in water vapor, and substantial dust settling 
to the disk midplanc appears to have already taken place. 
An initial estimate of the ortho/para ratio of the water 
suggests a value consistent with a high temperature equi- 
librium (0/P=3), but higher signal-to- noise observations 
of more lines is required to confirm a high O/P ratio. If 
confirmed, this is indicative of water formed via hot gas- 
phase chemistry, as opposed to water initially formed as 
ice in the outer disk and transported inwards as part of 
migrating solids. 

The spectra demonstrate that it is possible to obtain 
high resolution spectra of individual rotational lines of 
water in protoplanetary disks from the ground, in spite 
of telluric water absorption. These new tracers of the 
warm molecular layer have several unique advantages: 
The excitation energies and collisional rates of the high 
rotational levels are very similar to those of the ubiqui- 
tous fundamental ro-vibrational band of CO at 4.7/im. 
Given the chemical stability of CO, this allows for de- 
tailed studies of the relative water abundances in a well- 
defined location in the disk surface. The lines are strong, 
also relative to the continuum, and their analysis is not 
hampered by strong photospheric lines from the central 
star - in contrast to the weaker, much higher energy wa- 
ter lines in the 3 /xm "hot band" (Salyk et al. 2008)! The 
A-band lines may provide our only opportunity to ob- 
tain detailed velocity information of water at — 1 — 10 AU 
in disks around young solar analogs for the foreseeable 
future. Their observability from the ground also make 
them available to the very high spatial resolution line 
imaging of 10m class telescopes, as well as the next gen- 
eration of Extremely Large Telescopes. As such, they 
will be important complements to the lower resolution 
spectroscopy of the James Webb Space Telescope. 



Based on observations made at the ESO Paranal Ob- 
servatory under program IDs 084.C-0635 and 179. C- 
0151. The authors acknowledge valuable discussions with 
Alain Smette. 
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